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Abstract 

We investigate effects of the application of a kick to one-dimensional matter- wave solitons in a self-attractive Bose- 
Einstein condensate trapped in a optical lattice. The resulting soliton's dynamics is studied within the framework 
of the time-dependent nonpolynomial Schrodinger equation. The crossover from the pinning to quasi-free motion 
crucially depends on the size of the kick, strength of the self-attraction, and parameters of the optical lattice. 
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1. Introduction 

It is well known that Bose-Einstein condensates 
(BECs) with repulsive inter-atomic interactions, 
characterized by a positive scattering length, as > 
0, can give rise to stable localized matter-wave 
states in the form of gap solitons, in the presence 
of an optical lattice (OL), which induces a periodic 
potential acting on atoms. The existence of gap 
solitons in BEC was predicted theoretically [1] and 
demonstrated experimentally [2] , see also review [3] . 
Gap solitons are represented by stationary solutions 
to the respective Gross-Pitaevskii equation (GPE), 
with the eigenvalue (chemical potential) located in 
a finite bandgap of the OL-induced spectrum. Note 
that a periodic potential emulating the OL can also 
be induced by the spatially-periodic modulation of 
the transverse trap which confines the condensate 
in the transverse directions |4l5j . 

In the limit of a very deep OL, the underlying GPE 
can be reduced to an effective ID discrete nonlinear 



Schrodinger equation with the repulsive on-site non- 
linearity, assuming that the condensate's wave func- 
tion is a superposition of functions localized at par- 
ticular lattice sites, with a nearly vanishing overlap 
between them "^ . In this limit, gap solitons go over 
into staggered discrete solitons, which feature the al- 
ternation of the sign of the wave field between adja- 
cent sites of the lattice [7 . On the other hand, if the 
corresponding OL is weak, the general wave function 
may be split into a superposition of right- and left- 
traveling modes, giving rise to a system of coupled- 
mode equations for the respective slowly varying 
amplitudes (see, e.g., Ref. [8j), which is tantamount 
to the well-known coupled-mode equations for op- 
tical waves propagating through the Bragg grating 
[^. Actually, gap solitons were first predicted [10] 
and experimentally created [llj as optical pulses in 
fiber Bragg gratings. 

While the gap solitons are typically studied as 
standing localized states, a detailed numerical anal- 
ysis has demonstrated that an initial kick (repre- 
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sented by a linear phase profile imprinted into the 
soliton) may set that them in persistent motion, if 
the norm of the soliton (the number of atoms bound 
in it) is below a certain threshold value. The exis- 
tence of moving gap solitons was predicted in both 
one- and two-dimensional (ID and 2D) settings ^8,. 

In the BEC with attractive interactions (as < 0), 
solitons realize the ground state of the condensate. 
Such solitons were created in condensates of ^Li [12] 
and *^Rb [13j atoms, with the sign of the atomic 
interactions switched to attraction by means of the 
Feshbach-resonance technique (in the latter case, 
the solitons were observed in a post-collapse state of 
the condensate). In the presence of a periodic poten- 
tial, such solitons should exist too, with the chem- 
ical potential falling in the semi-infinite gap of the 
spectrum, as first shown in the context of the opti- 
cal setting [H] , and later demonstrated in detail in 
the framework of GPEs |15I16I5J . 

The objective of this work is to predict possi- 
ble regimes of motion of the ID matter-wave soli- 
tons in the self-attractive condensate, in the pres- 
ence of the OL potential. This will be done in the 
framework of the nonpolynomial Schrodinger equa- 
tion (NPSE), which is known to provide for high ac- 
curacy in the description of dynamical properties of 
solitons supported by the self-attraction. This model 
is described in Section II, while Section III summa- 
rizes the findings. The solitons will be first found as 
stationary solutions to the NPSE, and then set in 
motion by suddenly kicking them. Three dynami- 
cal regimes will be identified by systematic simula- 
tions: steady motion, gradual decay of the moving 
solitons, and firm pinning, when the kick cannot es- 
sentially disturb the initial soliton. In particular, re- 
gions where these regimes occur will be charted in 
a parameter plane of the model. The paper is con- 
cluded by a brief summary in Section IV. 

2. Bright solitons in the optical lattice 

Dynamics of atomic matter waves in rarefied BEC 
is very accurately described by the time-dependent 
three-dimensional GPE, 
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where ^(r,t) is the macroscopic wave function of 
the condensate (normalized to the total number of 
atoms, N), m the atomic mass, and Os the s-wave 
scattering length of the inter-atomic potential. As 



the trapping potential, we adopt the usual combina- 
tion of the tight confinement in the radial direction 
and longitudinal OL, 



U{r) = {l/2)muj'i{x^ + y^) + V{z), 
V{z) ^-Vq cos {2kLz), 



(2) 
(3) 



with /ci = 27r/A, where A is the wavelength of 
counterpropagating laser beams whose interference 
creates the OL, and Vb is its effective depth. In 
the subsequent analysis, we use the effective one- 
dimensional NPSE, which can be derived from Eq. 
((!]) by averaging in the transverse plane [29 . The 
NPSE has been demonstrated to be very accurate 
in reproducing results that can be obtained from 
the full 3D GPE; in particular, this approach has 
been tested for bright 17] and dark solitons [18], 
two-component condensates [20], and for the con- 
densate in a toroidal trap [TH], as well as for states 
with axial vorticity [21" . A similar approach, which 
generalizes the Thomas-Fermi approximation and 
may give still more accurate results, but does not 
apply to solitons, was later developed for tightly 
confined self- repulsive condensates |22j . 

The NPSE including axial OL potential ^ was 
derived too |16|5] , assuming the following factoriza- 
tion of the 3D wave function in Eq. ([1]), 
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where ID wave function /(z, t) is subject to normal- 
ization 



\I{z)\^dz 



(5) 



the adimensional interaction strength is g = 
2|as|iV/a_L, with a±_ = y^h/{muj±) the transverse 
trapping size. The respective NPSE takes the form 
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and 



with length and energy measured in units of a j 
?iaj_L, respectively. 

If applied to the ring-like (toroidal) geometry, Eq. 
© predicts that there exists a critical value, gmin, 
above which the axially uniform state in the torus 
becomes modulationally unstable, and is replaced, 
as the ground state, by a soliton-like configuration. 



Moreover, there is another critical value of the inter- 
action strength, gcoii, above which the localized so- 
lution ceases to exist due to the collapse, hence the 
soliton persists in interval gmin < 5 < 5coii- Similar 
effects in the toroidal setting were earlier studied in 
Refs. (23]. 

Both the full three-dimensional GPE and its 
NPSE reduction are mean-field equations. Contrary 
to the repulsive case, in the attractive case that 
we are dealing with, there is no superfluid-Mott 
insulator transition. Beyond-mean-field effects are 
irrelevant in our case too, since, as shown in Ref. 
[24], they appear at lattice sites populated with a 
small number of particles (< 5). In our case, due to 
the attractive interaction, most of the atoms reside 
in few highly populated sites. 

As concerns quantum effects, it may be relevant 
to mention that, unlike regular solitons in attractive 
condensates, gap solitons in the BEC with repulsion 
cannot represent the ground state. Therefore, in the 
framework of the full quantum description, the sta- 
bility of such states needs further analysis. In partic- 
ular, it was demonstrated that dark solitons in BEC 
with repulsive interactions between atoms are sub- 
ject to quantum depletion, with incoherent atoms 
gradually filling the dark-soliton's notch under the 
action of anomalous fluctuations [25,. Nevertheless, 
dark solitons, which accurately obey the mean- field 
description provided by the GPE, were successfully 
observed in the experiment [26]. Because loosely 
bound gap solitons may feature a set of notches in 
their "tails" , they may also be subject to the quan- 
tum depletion, in that sense. In terms of the quan- 
tum theory, this remains an open question, but the 
experimental results [213] clearly demonstrate that 
gap solitons can be created, following the predic- 
tions of the mean-field theory. On the other hand, 
for tightly bound soliton in the attractive BEC (as 
well as for tightly gap solitons in the repulsive con- 
densate, if they are created not too close to edges 
of the respective bandgap 127] ) , the quantum deple- 
tion should not be an issue, as their existence is not 
predicated on the presence of notches. 



3. Phase imprinting and ensuing dynamics 

As said above, the objective of this work is to an- 
alyze the dynamics of moving bright solitons in the 
model introduced above, which includes the peri- 
odic OL potential. The motion will be initiated by 
sudden application to the soliton of a kick with mo- 



mentum p. 

First, we generate a stationary soliton by numer- 
ically solving the corresponding stationary equa- 
tion, obtained from Eq. ([6]) by the substitution of 
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(z), where fx is the real chemical 



potential and (l){z) a real stationary wave function: 
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We fix fci = 1 in potential ([3]), and solve Eq. ([7]) 
with periodic boundary conditions, f{z) = f{z -\- 
L), with the total length of the ring-shaped system 
equivalent to 32 periods of the OL, L — 32tt. Tak- 
ing Vo = 0.5 for the OL depth, we find the ground 
state of Eq. ([7]) by means of the numerical scheme 
described in Refs. |16|5j . for values of the interac- 
tion strength g belonging to the above-mentioned 
existence range, ffmin < 9 < ffcoii- 

To initiate the dynamics, we multiply the so de- 
termined stationary solution (j){z) by exp(ipz), i.e., 
use the following initial configuration. 



/(z, t = 0) = 0(z) exp(ipz) , 



(8) 



where momentum p is compatible with the peri- 
odic boundary conditions, taking values p — n/16 
with integer n. We pick up such values in interval 
< p < 0.5. Configuration ([8]) can be created in 
the experiment by means of the well-known phase- 
imprinting technique |30j . The full time-dependent 
NPSE, Eq. ([6]), with initial condition ([U), was solved 
using the Crank-Nicholson predictor-corrector algo- 
rithm in real time [31j . 

Examples of the simulated evolution are displayed 
in Fig. [U which suggests the existence of different 
regimes corresponding to different values of interac- 
tion strength g. In particular, for small g (g = 0.15), 
the stationary soliton extends over many lattice cells 
and, after receiving the kick, it moves in quite a reg- 
ular way. At intermediate values of the interaction 
strength, g — 0.5, the soliton occupies a few sites, 
and its dynamics generated by the action of the kick 
is irregular. For a larger strength, g = 0.8, the initial 
wave packet is essentially localized in a single cell of 
the lattice, and stays trapped in the same cell (for 
the entire period of our simulation) after the appli- 
cation of the kick. 

To characterize the motion of the kicked soliton in 
a more accurate form, we calculate the average axial 
position of the soliton, zo(t) — J dz z |/(z, t)\ , and 
its average squared width. 
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Fig. 1. Profiles |/(^,i)| observed in the course of the evolu- 
tion of the kicked soliton, with initial momentum p = 0.25. 
From top to bottom: g = 0.15, g = 0.5 and g = 0.8. Here 
and in other figures, parameters of the optical lattice are 
k]^ = 1 and (unless specified otherwise) Vq = 0.5. 
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where the integration is performed over spatial pe- 
riod L = 327r (the expressions are not divided by 
the norm of the wave function because it is fixed 
to be 1, as per Eq. (O). These characteristics are 
shown, as functions of time, in Fig. [2] for the same 
values of g as in Fig. [1] The figure shows that, with 
g = 0.15, the soliton's center of mass (zq) moves at 
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Fig. 2. The center of mass, 20, and average axial width, 
(z^)^", of the soliton as functions of time. The initial mo- 
mentum is p = 0.25. 

a constant speed. We have verified that the respec- 
tive speed is well fitted by expression v = p/m*^ 
with the effective inverse mass of the ground state, 
1/m* = 1 — Vq /2, as predicted by the usual theory 
of the quasi- momentum in periodic lattices [32j . Si- 
multaneously, the width of the soliton, ^{z'^), dis- 
plays periodic small-amplitude oscillations, i.e., the 
moving soliton behaves as a breather. At g = 0.5, the 
motion of the soliton is less regular: its velocity is 
not constant, and the width increases in the course 
of the evolution, i.e., the soliton gradually spreads 
out; however, after becoming much broader than it 
was initially (by a factor ^^ 6), it seems to stabilize 
itself. At g = 0.8, the initial width of the soliton, 
^y{z^t = 0)), is, as said above, comparable to the 
size of the lattice cell. In this case, the soliton's cen- 
ter of mass does not exhibit any systematic motion, 
while the width initially increases (by a factor ^ 4) , 
but then ceases to grow. 

Thus, we do not observe complete delocalization 
in the regime of spreading . Instead, it is concluded 
that the soliton interacts with the OL in a complex 
way, emitting small-amplitude waves, which affect 
the time evolution of the average width. Figure [3] 
shows the soliton's profile a,t t = 108 in the case 
corresponding to the central panel of Fig. [U which 
features the maximum value of ■^/(z^). The figure 
demonstrates the presence of "splinters" , to the left 
of the soliton, which actually break off from the soli- 
ton shortly after the application of the initial kick. 
They travel backwards with respect to the soliton, 
reaching the largest distance from it at i = 108, due 
to the periodic boundary conditions. In this situa- 
tion, expression ^ does not provide an adequate 
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Fig. 4. The same cases as in Fig. [2] arc shown here in terms 
of the evolution of entropy mOI I. 



Fig. 3. Snapshot of \f{z,t)\ at t = 108 for the soliton in 
the central panel of Fig. [T] which corresponds to the largest 
value of {z^)^!"^ 

measure to estimate the spreading of the sohton. 
In principle, after performing a round trip in the 
toroidal trap, the splinters may hit the soliton, but 
the study of this issue require extremely long simu- 
lations, which we did not perform in the framework 
of the present work. 

As a better measure of the spreading, we introduce 
the following effective entropy, 



5(t) = -5]^„(t) 



In(A^coii) 



(10) 



where A„(i) is the share of the norm located, at 
time t, within the n-th lattice cell, i.e. An{t) = 
/Jn dz\f{z, t) 1^, and A^ccii is the total number of 
cells [recall we run the simulations for A^ccii = 32, 
and X]n=i' ^n ~ ^ ' ^^"^ ^° normalization condition 
(O]. Definition pil)) yields the maximum of the en- 
tropy, S — 1, if the matter is distributed uniformly, 
f{z,t) — const, and the entropy attains its mini- 
mum, S* = 0, if the entire norm is concentrated in a 
single cell. 

We have found that the unnormalized effective 
entropy, i.e. ^(t) In(A^ccii), turns out to be practi- 
cally independent of the total length of the system 
in the axial direction. In particular, while the mean- 
square width is affected by the backscattered ra- 
diation emitted after the kick, the entropy is not. 
More precisely, we have checked that the result of 
the computation of the unnormalized entropy does 
not change with the variation of the total lengths. 

In Fig. [4] we plot S{t) as a function of time for 
the same runs which were included in Fig. [2l It is 
seen that, for 5 = 0.15, initial entropy 5(0) is large, 
because the soliton covers many cells, and, in the 
course of the evolution, S{t) exhibits periodic small- 
amplitude oscillations around S{0). On the contrary. 



for g = 0.5 the stationary soliton covers few cells, 
hence initial entropy S{0) is quite small, but the 
entropy increases as a consequence of the soliton's 
spreading out. For g = 0.8, the initial entropy is 
similar to that in the case of g — 0.5, but its evo- 
lution is completely different, featuring very small 
aperiodic oscillations around 5(0), i.e. S{t) remains 
nearly constant, as the narrow strongly pinned soli- 
ton does not start conspicuous motion or spreading 
out after being kicked. 

Thus, our systematic simulations reveal the exis- 
tence of three different dynamical regimes: (i) sta- 
ble breathers, i.e., solitons steadily moving across 
the lattice potential at an almost constant velocity, 
with small-amplitude shape oscillations; (ii) disper- 
sive dynamics, in which case the soliton strongly 
spreads out in the course of the evolution; (iii) local- 
ization, in which a narrow soliton remains trapped 
in one lattice cell. In Fig.[5l we plot the respective di- 
agram in the parameter plane of (p, g) , where these 
three regimes are mapped as follows. 

(i) The black region represents steadily traveling 
solitons, whose effective entropy varies by < 5% in 
the course of the long evolution. 

(ii) The gray region: the dispersive regime, in 
which the solitons move at a variable speed, and 
their effective entropy increases by more than 5% 
against the initial value. 

(iii) The white region: localized solitons, with the 
center of mass firmly pinned to the initial position. 
In the latter case, the solitons may be slightly dis- 
persive at the initial stage of the evolution, but their 
effective entropy is always much smaller than in the 
other two cases. 

Figure [5] summarizes results of 560 numerical 
runs, with 16 values of p and 35 values of g. The 
momentum resolution is imposed by the periodic 
boundary conditions: Ap = 2Tr/L, where L = 327r 
is the total axial length, as defined above. By choos- 
ing larger L, one can reduce Ap. The step in the 
variation of strength g in Fig.[5]is Ag — 0.25. 
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Fig. 5. The map of different dynamical regimes in tlie plane 
of {p,g), where p is the momentum imparted to the soliton 
by the initial kick, and g the strength of the self-attractive 
nonlinearity. The black region: steadily moving breather-like 
solitons; the gray region: spreading out of irregularly moving 
solitons; the white region: firm pinning (the soliton does not 
start progressive motion). Lattice parameters: Vq = 0.5 and 
kL = 1. 



Figure [H] shows that the regime of the stable mo- 
tion is confined to the region of small values of inter- 
action strength g. For fixed g, a transition from the 
steady motion to the dispersive regime is observed 
with the increase of kick p, and the numerical results 
suggest that the steady-motion region disappears as 
p attains large values. As expected, the increase of 
g drives the system into the regime of pinning, but 
one can leave it by increasing the kick. 

It should be stressed that, while the border be- 
tween the localization and other two regimes is well 
defined, the exact location of the boundary between 
the steady-motion and dispersive regimes depends 
on our choice of the 5% maximum for the allowed 
change of the entropy in the course of the evolu- 
tion. Admitting a 10% change does not conspicu- 
ously affects the location of the border, but defining 
the threshold of 20% makes the region of the stable 
motion slightly larger (not shown here in detail) . 

At the boundary between motion and dispersion 
regimes, we observe a density profile that becomes 
more and more noisy, while the velocity of the trav- 
eling soliton decreases in an irregular way, similar to 
some regimes of the motion of gap solitons in the ID 
model with the repulsive cubic nonlinearity, which 
were reported in Rcf. '8^ . This drop in the velocity is 
also relevant to the identification of the boundary of 
the pinning regime, where, typically, the kicked soli- 
ton visits a few lattice sites and then stops. For in- 
stance, at g — 0.55 and k — 0.25, the soliton passes 
a few sites, then bounces back to the original posi- 



Fig. 6. The same as in Fig. [S] but for Vo = 1. 

tion, and stays pinned there. 

The diagram of the dynamical regimes depends 
on parameters of the optical lattice, Vo and fci, see 
Eq. ([3]). In particular, Fig.[6l obtained by collecting 
results of other 16 x 33 = 528 numerical runs, dis- 
plays this diagram for a stronger lattice, with Vo = 
1, while the OL period remains equal to tt (i.e., fc^ = 
1). Qualitatively, the picture is similar to that shown 
in Fig. [5l but there are apparent quantitative dif- 
ferences. In particular, the regions of the stable mo- 
tion and dispersive behavior are strongly reduced, 
which is not surprising, as the deeper OL can pin the 
solitons stronger. In addition, the figure shows that 
there exists a critical momentum, pc, above which 
the region of steady motion does not exist anymore. 
This critical value pc decreases with the increase of 
OL depth V"o. 

As mentioned above, the upper bound for the con- 
sidered values of the interaction strength, g, is im- 
posed by the onset of the collapse, g = gcoU- In par- 
ticular, in Ref. [16 it was found that gc = 1.07 for 
Vo = 0.5, and g^ou = 0.96 for Vq = 1 (if /cl = 1 
is fixed). On the other hand, the minimum strength 
above which soliton states exist is gmin ~ tt^/L [23] : 
in this work, this border does not manifest itself, as 
we considered values of g at which the soliton is well 
localized within the lattice; in fact, we always took 
./mO))<L/lO. 

In an attempt to understand the behavior of the 
kicked soliton in analytical terms, we tried a vari- 
ational approach. Using a Gaussian ansatz, we de- 
rived the respective variational equations of motion. 
We have thus obtained a set of equations equiva- 
lent to that recently reported in Ref. [33] ■ However, 
a careful analysis demonstrates that this approach 
fails to describe the soliton's dynamics correctly in 
the present model: although it does predict regions 



of pinning, moving breathers, and solutions whose 
width expands indefinitely, the resulting diagram is 
quite different from those in Figs.[5]and[6l This prob- 
lem might be expected, since, even in the static case, 
the Gaussian ansatz does not reproduce the behav- 
ior revealed by numerical solutions ^16j . 

4. Conclusions 

We have studied the dynamics of moving matter- 
wave solitons in the self-attractive BEG trapped 
in the axial optical-lattice potential. The dynamics 
was induced by kicking a stationary soliton via the 
phase imprinting. We have performed systematic 
simulations of the time-dependent nonpolynomial 
Schrodinger equation, which accurately describes 
Bose-Einstein condensates under the tight trans- 
verse confinement. 

Three dynamical regimes have been identified. 
The first of them is the steady motion of stable 
solitons with small-amplitude intrinsic oscillations, 
which is observed at relatively small values of the 
strength of the inter- atomic attraction g, and rel- 
atively small size of the kick p, and the second is 
the dispersive regime, in which the soliton's veloc- 
ity decreases irregularly, while the soliton suffers 
systematic spreading. The latter outcome of the 
application of the kick to the soliton also occurs for 
relatively small g, but at larger p. The third regime 
naturally features firm pinning, at large values of 
g and small p. A somewhat similar description for 
dynamical regimes was elaborated in Ref. [6], but 
the position of the different regions in the respec- 
tive parameter plane is totally different, due to 
the fact that the quasi-discrete limit considered in 
that work is not adequate for comparison with the 
present model. 

Available experimental techniques [12113] make 
the experimental realization of the predicted dy- 
namical regimes quite feasible. It may also be rel- 
evant, for the purposes of the theory and plausible 
experiment alike, to study in details collisions be- 
tween steadily moving solitons. The latter problem 
will be considered elsewhere. 
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